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Abstract 

In this short note, we compute, for large n the determinant of a class 
of n x n Hankel matrices, which arise from a smooth perturbation 
of the Jacobi weight. For this purpose, we employ the same idea 
used in previous papers, where the unknown determinant, D n [w a) ph\ 
is compared with the known determinant D n [w a ^]. Here w a> p is the 
Jacobi weight and w a ^h, where h = h(x), x S [—1,1] is strictly 
positive and real analytic, is the smooth perturbation on the Jacobi 
weight w a ^(x) := (1 — x) a (l + x)@ . Applying a previously known 
formula on the distribution function of linear statistics, we compute 
the large n asymptotics of D n [w a ^h] and supply a missing constant 
of the expansion. 
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1 Introduction and Preliminaries. 



The purpose of this note is to find heuristically an asymptotic expansion for 
determinants of certain Hankel matrices. The matrices are generated by the 
moments of a function defined on the interval [—1,1]. Let w(x) be a function 
of the form 

Wa,p( x )H x ) 

where 

w a>p (x) = (1 -x) a (l + xf, a>0,P>0 

and h(x) is a strictly positive function with a derivative satisfying a Lipschitz 
condition. 



Define 



/•ifcM = J x k w(x)dx, k = 0,1,2,.. 

and 

D n [w] = det^+kH)^. 

The motivation for investigating such perturbed Hankel determinants, comes 
from Random Matrix Theory and its applications, where one studies the 
generating functions of linear statistics |U |3] . Also see [I] and some of the 
references in that volume. 

Our goal will be to show formally that, for w = w a ^h, 



(1.1) D n [w] ~ 2-^+" + %^ 2 )/ 2 - 1 / 4 (27rrexp ( - C ^M= dx ) C 



7T J_ x VT 



X 



2 



where the n independent constant C is given by 



exp 



1 lnft w l P f^EZm iy]ix 



4vr 2 J_ x Vl - x 2 \ J -i y - x h(y) 



a+l f 1 \nh(x) \ G 2 (s±f±l) G 2 + l) T (gfg) 
X( ' X|>! 2tt J_ 1 ^/T^2 X ) G{a + /3 + l)G{a + l)G(J3 + l) 



2 



In the above formula the function G is the Barnes G-function, an entire 
function that satisfies the difference equation G(z + 1) = T(z)G(z), with 
G(l) = 1. The result (1.1) is also valid for a > —1/2, and (3 > —1/2, since 
this expression is real analytic in a and (3. 

The main idea, which can be traced back to a paper of Szego |3j, is that 
one can find the above formula in two steps. The first is to consider the 
"pure" weight 



Using some basic results from the theory of orthogonal polynomials the Han- 
kel determinant for the pure weight can be found exactly and then easily 
computed asymptotically. This step is rigorous and in fact may be the first 
instance where these asymptotics are found completely. 

The next step is to use the linear statistics formula derived from the 
Coulomb fluid approach [U E] expected to be valid for sufficiently large Ti- 
to compute the quotient 



thus achieving the desired result. 

We note that in a recent work [S], the asymptotic formula for D n appears, 
but without the constant term. In future work, we hope to use the techniques 
of |2j to make the ideas presented here complete and thus firmly establish 
the validity of the asymptotic formula. 

We begin with some notation. Let P n (x) be monic polynomials of degree 
n in x and orthogonal, with respect to a weight, w(x), x G [a, b]; 



w a>p (x) = (l-x) a (l + xf. 



(1.2) 



(1.3) 




where hj [w] is the square of the L 2 
respect to w, over [—1,1]. 



norm of the polynomials orthogonal with 
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From the orthogonality condition there follows the recurrence relation, 
(1.4) zP n (z) = P n+ i(z) + a n P n (z) + PnP n -i(z), n = 0, 1, 

where (3qP-i{z) := 0, a n , n = 0, 1, 2, ... is real and /3 n > 0, n — 1, 2, ... 

There is an intimate relationship between the values of /3 n , h n and the 
Hankel determinants. 

Indeed, the determinant, for any weight w, 

n-l 

D n [w] = U h AM- 

3=0 

In addition, 

3 

hj[w] = h [w}Y[(3 k . 

k=l 

Thus if we can compute /3j it follows that both hj and D n can be explicitly 
determined. For this and all other basic results see JU] • 

For the monic Jacobi polynomials, that is in the case when w = w a ,p, it 
is well-known that 

(3 2 -a 2 

Otr. 



n (2n + a + (3)(2n + a + (3 + 2) 
and 

4n{n + a){n + (3){n + a + (3) 

Pn 



{2n + a + P) 2 (2n + a + (3 + l)(2n + a + /3 - 1) ' 
Hence it follows that 



1.5) K[w a> p] 



2 2n +a+ /3+i r(n + l)T(n + a + l)T(n + (3 + l)T(n + a + f3+l) 
(2n + a + (3 + l)[r(2n + a + (3 + l)] 2 
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and 



1.6) D n [w a> p\ 

r ( a+(3+l \ p2 (a+£+l\ pi ( a+j3 , •. \ 
_ r>-n(n+a+l3)fr> \n L \ 2 ) _ V 2 I _ I 2 1 ^ 

1 J Gf(a + /? + l)G(ar+l)G(/?+l) 
G(n + l)G(n + a + l)G(n + /3 + l)G(n + a + /3 + 1) 



X G 2 („ + a+l+1) G 2 („ + «±g + ^ r + 

where G(z) is the Barnes G-function. See jH] for a first-principle derivation 
of the recurrence coefficients. 

The asymptotics of the Gamma function and the Barnes G— function are 
well understood. We have that 



r(n + a) ~ v / 2^e-"n n+a - 1/2 ) 

G (n + a+l) ~ n (n + a)V2-l/12 e -3nV4-an (27r) (n+a)/2^ 

where K := G 2 / 3 (1/2)tt 1 / 6 2- 1 / 36 . 
^From the above asymptotic expressions an easy computation shows that, 

D n [w a ,p} ~ 2 -"( n+Q+ ' 3 )n (Q2+/32)/2 - 1/4 (27r) n 

G 2 (o+i+i) G 2 («±e + 1) r (^±f+i) 



;i-7) 



G(a + /3 + + + l; 



The above formula is the promised result for the "pure" weight. When a 
= P, we find, 



7T n I 

" 1 -y2 / 



A*Ko] ~ ^i^rry G 2 (l/2)r(l/2). 

This is consistent with Hilbert's [Zj asymptotic expression for large n, of the 
Hankel determinant associated with the Legendre weight, 

(D n [w ,o}) 1/n = 7 —— (1 + e n ) , where lim e n = 0. 

Z n—>-oo 

and we have changed the notations of 7\ to be compatible with ours. 
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2 Perturbed Jacobi Weight 



In this section we show how to compare the unknown Hankel determinant 
D n [w at ph], with the known Hankel determinant D n [w a ^]. It is known from 
PU| (see also jHE]) that 

(2.1) ^4=/f[^) 

where = Q(xi, x n ), and 

^ , i _ A - A Q Ui< J<k <n( x k ~ xj) 2 Iir=i w^jx^dxj 

A -An 

This can be rewritten as an average of the exponential of the linear statis- 
tics Ya=i ^K x i), i- e -> 

( 7 

exp j In h(xi) 

Note, because of the assumptions on h, Inh is well defined for x G [—1, 1]. 
Results, at least in a heuristic way, are known about such linear statistics. 
In particular, the logarithm of (2.1) is, for large n, 



1 A 1*M*) (p [ bn V(bn-y)(y-a n )h>(y) d L 



4ff2 A y/(b n -x)(x-a n ) \ J an V-x h(y) 

+ / \nh{x)a(x)dx, (2.3) 

J a„ 

where the equilibrium density a(x), defined for x G [a n , b n ] is 
a/(&„ - a;) (a; - a n ) f bn v'(x) - v'(y) dy 

(Tlx) = — ' 



27r2 A x-V \/{b n -y){y-a n y 
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and 

v [x) ■— 



yj a ,p(x) x — l x + 1 

The end points a n , b n , of the support are determined by, 

: K xv'(x) 

2im = / — ax, 



n s/Wn -x)(x - a n ) 
= / v ; dx. 

Ja n V(6 n - X){x ~ On) 

The equation (2.3) was derived in [3] and later a large n version was found 
in [Sj. For another derivation of (2.3) using a "small fluctuations" approach 
see [IJ. 

In our problem, the above equations become, 

a + 3\ a 3 

n+\— r— = — = = + 



2 J 2y/(l-a n )(l-b n ) 2y/(l + a n )(l+b n ) 

= ° 2 

y/(l-a n )(l-b n ) y/(l + a n )(l+b n y 

and the solutions are 



3 2 - a 2 - 4A/n(n + a)(ra + 3){n + a + 3) 
(2n + a + 3 + 2) 2 



/? 2 - a 2 + 4^72(72 + a)(ra + /3)(n + a + 3) 
n ~ (2n + a + /3 + 2) 2 ' 

In the Coulomb fluid approximations the diagonal (a„) and off-diagonal 
recurrence coefficients ($ n ) are 

b n + a n (3 2 - a 2 



A, 



2 " (2n + a + /3) 2 ' 
(b n — a n ) 2 4n(n + a)(n + 0){n + a + 0) 
16 ~ (2ra + a + /?) 4 ' 
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and the deviations from the exact results are, 



Pn-Pn = 



16n 2 



1 +0 (X\. 



rr 



For later reference we also note that, 



1 + a n = ^+Ol — 



2n 2 



TV 



' + o' ' 



2n 2 



A simple calculation shows that, for x G [a n , fe n ], 



aix) 



\J(p n -x){x- a n ) 



1 

2^ 
= | n + 



P 



^(l-a n )(l-b n )(l-x) v /(i + an )(i + 6 n )(i + 



.t; 



a + /9 



1 — x- 



-, -1 < a n < fo n < 1. 



where we have used, 



a 



V(l-On)(l-6n) 



n + 



= n + 



a + /3 
a + /3 



v/(l + a n )(l + 6„) 2 
Therefore, for x G (—1, 1), and n large, 



-(-) " + + 
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Put f(x) = \nh(x), and x = R n + r n t, where R n := (b n + a n )/2, and r n :- 
(b n — a n )/2, the second term of (2.3) becomes 



a + (3 

n + 



2 f 1 f(R n + r n t) 



2 

while the first term of (2.3) reads, 

r„ f 1 f(Rn + r„s) ( [> vTHlfVn 1 , 

Now, since 

K = ^i + °(4 

4 7Jr \n d 



we see that the second term of (2.3) is asymptotic to 



a 



+ p\ f 1 lnh(x 



n-\ 7r —) I — v ' = dx + o(l) 



2 J J_ x nVl 



x 



while the first term of (2.3) is asymptotic to 



1 inM*) ( P rvnzm iy u x+om . 



4vr 2 J_ x y/l - x 2 y y_i 2/ - x %) 
The above two expressions combined with (1.7) give the formula (1.1). 
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